Rotating Bose-Einstein condensates with attractive interactions 
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We examine the phase diagram of a Bose-Einstein condensate of atoms, interacting with an 
attractive pseudopotential, in a quadratic-plus-quartic potential trap rotating at a given rate. In- 
vestigating the behavior of the gas as a function of interaction strength and rotational frequency of 
the trap, we find that the phase diagram has three distinct phases, one with vortex excitation, one 
with center of mass excitation, and an unstable phase in which the gas collapses. 
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I. INTRODUCTION 

The velocity field of a superfluid is irrotational. As 
a result, such fluids develop vortex states when rotated, 
and the circulation of the superfluid velocity around any 
closed path is quantized Q . Trapped Bose-Einstein con- 
densates of alkali-metal atoms provide ideal systems for 
testing these ideas. Remarkably, they permit experimen- 
talists not only to vary the form of the trapping potential 
but also to change the sign of the scattering length asso- 
ciated with the effective atom-atom interaction. 

In a harmonic trap, vortices in Bose-Einstein con- 
densates with repulsive interactions are always singly- 
quantized UQ- If, however, the effective interaction is 
attractive or the trapping potential has a different shape 
, the situation can change dramatically. The recent ex- 
periment of Rcf. 4] investigated the behavior under ro- 
tation of an effectively repulsive Bose-Einstein conden- 
sate confined in an anharmonic trap. As shown in the 
theoretical studies of Refs. jSIESIslliEIIlIll, vor- 
tices can be multiply quantized in any trap that grows 
faster than quadratically. Further, Refs. [IsllJjUJl have 
demonstrated that the angular momentum is carried by 
center-of-mass motion in an effectively-attractive conden- 
sate that rotates in a harmonic trap. 

In this study we investigate the phase diagram of an 
anharmonically-trapped Bose gas in the case of an ef- 
fective attractive atomic interaction as a function of the 
rotation frequency of the trap and of the strength of the 
dimensionless coupling constant, Na/Z, where N is the 
total number of atoms, a is the s-wave scattering length, 
and Z is height of the gas along the rotation axis. Both 
the sign and the magnitude of the scattering length can 
be adjusted using Feshbach resonances. The magnitude 
of the coupling constant can also be controlled by chang- 
ing either the number of atoms in the trap or the trap fre- 
quencies. For sufficiently strong attractive interactions, 
such systems are unstable against collapse 0, 0, 
However, as long as the interaction is balanced by the 
zero-point motion of the atoms in the trap (and also by 
the kinetic energy associated with the vortices in the case 
of a rotating gas) , the gas can exist in a metastable state. 



In discussing the effects of rotation on the metasta- 
bility of a gas with attractive interactions we must dis- 
tinguish two physical situations, whether the trap is be- 
ing rotated at a given angular frequency - the principal 
case we consider here - or whether the gas is set into 
rotation, and then the rotation of the trap is switched 
off. The former case, the analog of the Meissner effect 
in superconductors, examines the equilibrium states of 
the system in the presence of rotation. The latter exam- 
ines the metastability of supercurrents in the absence of 
a driving force. As stressed by Leggett [l^, superfluid 
flow in a system with effectively attractive interactions is 
not stable in the absence of external rotation. Thus the 
rotating states that we find here to be metastable in a 
rotating trap should lose angular momentum and cease 
rotating once the rotation of the trap is switched off. The 
details of how the angular momentum associated with su- 
perfluid flow in such systems is lost is beyond the scope 
of this paper. Examples of the breakup of a vortex are 
given in Refs. [2(1 121| . Angular momentum can also be 
continuously taken out of the flow through interaction of 
the flow with normal fluid. 

Anharmonic traps allow a rich structure of three phases 
depending on the strength of the attractive interactions 
and on the rotational frequency. The first is an unsta- 
ble phase in which the gas collapses to a more dense 
state. The second phase is characterized by the angular 
momentum being carried by the motion of the center of 
mass. The third phase involves either a mixed state of 
multiply and singly quantized vortices or a pure state of 
multiple quantization. The general structure of the corre- 
sponding phase diagram is shown schematically in Fig. 1. 
References [22| have investigated this problem using both 
variational and exact numerical techniques. It is crucial 
that the trap be anharmonic to have such a structure. In 
a harmonic trap, repulsive interactions ensure that the 
critical frequency for rotation is smaller than the trap 
frequency and thus permit the formation of vortices Q|. 
Attractive interactions increase the critical frequency for 
rotation to a value larger than the trap frequency with 
the result that the centifugal force cannot be balanced by 
the restoring force of the trap and the atoms fly apart. 
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FIG. 1: Schematic phase diagram of a rotating Bose-Einstein 
condensate trapped in an anharmonic potential. The origin 
is the upper left corner. The vertical axis is the strength of 
the attractive effective atom-atom coupling constant and the 
horizontal axis is the rotational frequency of the trap. 



In analyzing the problem we consider atoms interacting 
via a short-range effective interaction, 



(1) 



Here, Uo = inh a/M is the strength of the effective two- 
body interaction, where M is the atomic mass. We con- 
sider an anharmonic trapping potential of the form 



V(p,z) 



1 



1 + A 



+ V z (z). (2) 



Here, u> is the trap frequency of the harmonic poten- 
tial perpendicular to the z axis (which is assumed to be 
the axis of rotation) and do = {Ti/Muj) 1 / 2 is the oscilla- 
tor length. The dimensionless constant A is small (e.g., 
A « 5 x 10 -3 in the experiment of Ref. Q) and V z is the 
trapping potential along the z axis. Denoting the excita- 
tion energy of the first excited state along the z axis by 
AE Z , we assume that 



nU < hu < AE Z 



(3) 



where n is the typical atom density. The inequality be- 
tween the left and the right terms implies that the cloud 
is in its lowest state of motion along the z axis, and the 
problem thus becomes effectively two dimensional. Fur- 
thermore, under the above conditions, the typical den- 
sity is ~ N/cLqZ. The interaction energy nllo is ~ aahuj, 
where a = N/Z is the atom density per unit length. The 
assumption of small nllo in the left inequality of Eq. J3J), 
is equivalent to ca < 1. 

In the following section, we examine the phase where 
the rotating gas forms vortices. In Sec. Ill we study 



the phase in which the angular momentum is carried by 
the center of mass, and in Sec. IV examine the unstable 
phase. Section V discusses the general features of the 
phase diagram, and Sec. VI summarizes our conclusions. 



II. VORTEX PHASE 

For sufficiently weak interactions, the behavior of the 
system under rotation is dominated by the anharmonic- 
ity of the trap, and the many-body wave function is a 
product of single-particle states describing multiply - or 
singly-quantized vortex states. Assuming weak interac- 
tion (a a <C 1) and weak anharmonicity (A <C 1), we can 
restrict our attention to a basis of eigenstates of the har- 
monic potential with zero radial excitations and angular 
momentum mh. (We now set h = M = w = 1 for con- 
venience. We denote many-body states by capital letters 
and single-particle states by small letters.) The basis 
states are 



1pm(P, <t>) 



1 



72 



(4) 



where m is a non- negative integer, and z = pe 1 ^ . The 
symmetrized mean-field many-body state with total an- 
gular momentum L can then be expanded in this basis 
as 



N oo 



* v (ri,...,rjv) = JJ^ 



(5) 



witn Sm=o TO l c ™l 2 = L/N, and J2m=o |c m | 2 = 1, where 
"Jo is the many-body state of the nonrotating cloud, 



*o = 



,N/2 



exp 



\Si\ 2 /2 (6) 



Here, 4> z is the ground-state many-body wavefunction of 
the cloud along the z axis. Henceforth, we neglect $ z , 
since it plays no role in our analysis. 

The energy of the system in the state \l/ v was calcu- 
lated in Refs. 0,^2- Here, we note only certain features 
related to the trapping potential. If w, = — R are rela- 
tive coordinates (with R the center of mass coordinate), 
the potential energy can be written as 



V(pi, ...,p N ) = 



JV 



£ K + m 4 ) 

i=l 

N 

NR 2 + N\R 4 + 4Ai? 2 £ 



(7) 



In the state ty v , the first terms, involving Wi, are of or- 
der NTiui and XNhto; by contrast the last three terms in 
Eq. J7J are of order Tiuj (i.e., unity), Xhw/N, and \hu>, 
respectively, since N(R 2 } is of order unity. Therefore, in 
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FIG. 2: The phase diagram of a Bose-Einstein condensate 
trapped in a quadratic-plus-quartic potential in the Q/lj - 
aa plane for A = 0.05. The vertical line denotes the critical 
frequency for exciting the center of mass (CM), and the lower 
dotted line gives the boundary between vortex and center- 
of-mass motion. To the left of the vertical line, the system 
does not rotate. The higher dashed line denotes the boundary 
between pure multiply-quantized vortex states with winding 
numbers m = 1 to 9 and a mixed phase consisting of multiply 
and singly quantized vortices. 



the frame rotating with angular frequency, Q, the energy 
per particle of a multiply-quantized vortex state, 



is 



K 

Nhuj 



*v(ri 



1 + m 1 



N 

n 



F*0, 



(8) 



A 



(m + l)(m - 



+aa- 



(2m)\ 
2 2m (m!) 



P>2 ; 



(9) 



as in Refs. [TSIH. 

As shown in these references, if cr\a\ exceeds some crit- 
ical value, the state containing a multiply quantized vor- 
tex, described by ip m , is unstable against becoming a 
state of mixed angular momentum of the form 



1p = C m -i4> m -l + Cmlpr, 



m+lV'm+l- 



(10) 



The dashed curves in Fig. 2 show this phase boundary 
with A = 0.05, for m = 1 to 9 from left to right. In each 
of the regions above the dashed lines, the energy of the 
gas in the rotating frame is minimized for \c m \ = 1, and 
\c m ±i\ = 0. 



III. CENTER OF MASS EXCITATION 

If the scattering length a is sufficiently negative, we 
expect the angular momentum to be carried by the center 
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FIG. 3: The same graph as in Fig. 2 for a wider range of values 
of aa. The higher curve separating the vortex phase and the 
center of mass phase is the lower curve shown in Fig. 2. 



of mass motion |0, 0] . The many-body wavefunction is 
then, 

*cm(ri,...,rjr) = (^^) *°' 

where L is non-negative. Calculating the energy per par- 
ticle of the system in the state f cm in the rotating frame, 
we find that 

§^ = l + l(l-^} + ^ +il + 2) + ( ra, (12) 

with I = L/N. The term A/ 2 /2 comes from the quartic 
component cx i? 4 in Eq. (JJJ, the term 2X1 results from 
the component X^ li wf, while the term A results from 
4:\R 2 J2i w i- The final term in Eg. ((T2"|l is the interac- 
tion energy of a non-rotating cloud. Since center of mass 
motion does not affect the relative coordinates, (Vint) is 
unaltered. 

Differentiating Eq. (|12fl with respect to I we find that 
the critical frequency of rotation of the cloud (for center 
of mass excitation) is [23 , 



tt c /uj = 1 + 2A, 



(13) 



which is shown, for A = 0.05, as the vertical curve in 
Figs. 2 and 3. 

For sufficiently negative interaction strength the en- 
ergy of f cm becomes smaller than that of ^v, and we 
expect a phase transition from the vortex phase investi- 
gated in Sec. II to the center of mass state. The overlap 
between these two states is exponential small in the num- 
ber of atoms. In the limit of a large number of atoms, 
this transition is discontinuous in contrast to the tran- 
sitions examined in Sec. II. In the case of a pure giant 
vortex given by Eq. JHJ, for example, the overlap is 



(*cm|*v) = 



N L (m\) N 



1/2 



(14) 



4 



with L — mN . In the limit N — > oo, 



(*cm|*v> = (2 7 rm) 1 / 2 
= (2 7 rm)- JV / 4 



N/2 



(15) 



where last equality applies for m>l. 

The phase boundary in the £1/lu - aa plane is approxi- 
mately the line along which E' v = E' cm , given by Eqs. © 
and (|12(l . This result is shown as the lower (higher) 
curve in Fig. 2 (Fig. 3). In fact, the true phase bound- 
ary lies even lower since the energy for the single-particle 
mixed states of the form ip = • • • c m _ii/) m _i + c m ip m + 
Cm+iipm+i • • ■ is necessarily smaller than that for ip = ip m 
in this region. 



IV. UNSTABLE PHASE 

Sufficiently large attractive interactions render the 
cloud unable to support itself, and a dense atomic state 
results. In a nonrotating system, this instability occurs 
when a\a\ is on the order of unity [Til Il7| . As we show 
below, in a rotating cloud this value is even more nega- 
tive because of the kinetic energy of the rotational motion 

M 

To calculate the phase boundary, we start with the 
energy in the rotating frame in the state f cm , Eq. (|12(l . 
but regard the corresponding oscillator length as a varia- 
tional parameter. We previously assumed that the oscil- 
lator length is fixed and given by {h/NMujf/ 2 . Here, we 
perform the same calculation assuming that the oscilla- 
tor length is ftifi/NMio) 1 / 2 , where (3 is real and positive. 
The result is 



Nhuj 



u 2~ • - [3 2 ' 

(16) 

For a non-rotating cloud, the value of /3 that minimizes 
the energy is 



1 + 2aa 
1 + 2A 



1/4 



(17) 



which implies that the critical value for collapse is aa 
Differentiating E' cm with respect to I we obtain 



I = 



2A/3 2 



1 

A/3 1 ' 



(18) 



so that the critical frequency for center of mass excitation 
is 



uj 2\ + (3 2 



2\(3 2 



(19) 



This expression reduces to Eq. (|13fl when (3 — > 1 in the 
limit of weak interactions, aa <C 1 and A <C 1. The point 



-0.2 



-0.4 



-0.6 



m = -1/ \ m = 


i 

m = 1 


\ Uniform / 




y Localized "--'' 





-1 -0.5 0.5 1 
Q 

FIG. 4: The phase diagram of an effectively attractive Bose- 
Einstein condensate confined in a toroidal trap, from Ref. |2^l . 
Here 7 is the ratio between the interaction energy and the 
kinetic energy, corresponding to aa in the present problem, 
and Q, is the dimensionless rotational frequency of the torus. 
The function ~y(Sl) is periodic. 



o o 




nr 



FIG. 5: Equipotential plots of the density (higher) and phase 
(lower) of the condensate in the state of the form of Eq. HUH 
for m = 3. In the left graphs C2 = C4 = 0, while in the 
right ones all c m are nonzero. The left graphs correspond to 
the uniform phase of the one-dimensional problem, while the 
right ones correspond to the localized phase. The axes extend 
from — 4do up to Ado- 



where E' cml expressed in terms of (3, Cl/u, aa, and A, 
ceases to have a local minimum as function of (3, signals 
the collapse of the cloud. The result for A = 0.05 and 
varying Vl/lo and aa is shown in the lower line of Fig. 3. 



V. FEATURES OF THE PHASE BOUNDARY 

The calculated phase diagram of Fig. 2 is universal 
[Tfl Il2l | in the sense that the same graph would be ob- 
tained for A' = aA by rescaling the axes by a, i.e., 
(1 — Q/uj) — > a(l — Cl/u) and aa — > aaa. Further, the 
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dashed line indicating the phase boundary between pure 
multiple quantization and the mixed phase is exact in 
the limit of small A and a\a | OHEl. In the phases of 
multiple quantization, the energy in the rotating frame is 
minimized (i.e., the energy is an extremum and its second 
derivative with respect to admixtures of other states is 
always positive). Thus, even though the effective interac- 
tion is attractive in this system, the vortices of the driven 
system are stable against small perturbations. Further- 
more, as one decreases Cl/ui crossing the segments of fixed 
m, there is a continuous path along which the energy de- 
creases towards the absolute minimum of the energy (in 
the lab frame) . This implies that there can be no persis- 
tent currents, in agreement with the toy model presented 
in Sec. VI of Ref. O. 

The higher phase boundary in Fig. 2 closely resembles 
that of a rotating Bose-Einstein condensate confined in 
a one-dimensional toroidal trap, with an effective attrac- 
tive interaction between the atoms. This system was in- 
vestigated in Refs. i its phase diagram is shown in 
Fig. 4. In this one-dimensional problem the instability 
of states of the form 4> m — e lm ^ / \/2tt is actually towards 
the combination c m _i0 m _i + c m (f> m + c m+ i<f> m+1 , which 
is the one-dimensional analogue of our problem. [See 
Eq. I|1U[1 .] The reason for this resemblance is the form of 
the effective potential felt by the atoms, 

V eS = V - MQ 2 p 2 /2, (20) 

which has a mexican-hat shape for Q > uj. We identify 
the phase of pure multiple quantization as the uniform 
state of the one-dimensional problem and the other two 
phases as the localized state of the one-dimensional case. 
References exclude the possibility of a state of the 
form of Eg. (|1(J|) . However, at least for the parameter 
range considered here close to the phase boundary, there 
is also a mixed phase for which the order parameter has 
the form of Eq. lfTU|) . 

When c m _i = c m+ i = 0, these states describe 
multiply-quantized vortices of winding number m, and 
the density is homogeneous, as shown on the left up- 
per graph of Fig. 5 (for m = 3.) We see in the lower 
left graph a triply quantized vortex at the origin. On 
the other hand, when c m _i and c m +\ are nonzero, these 
states describe a combination of a multiply-quantized 
vortex (m = 2) at the center of the cloud with wind- 
ing number m — 1, here 2 in the lower right figure, one 



singly-quantized vortex just to left close to the origin, 
and another not shown, much further from the center, 
where the density is very low. As shown in the upper 
right graph of Fig. 5 (for m = 3), the density is now in- 
homogeneous (and resembles that corresponding to the 
center-of-mass phase) . This is an energetically favourable 
configuration because the effective atom-atom attraction 
prefers the inhomogeneity around the minimum of the ef- 
fective potential V c s. The reason why the minima in the 
phase boundary of Fig. 2 decrease with increasing Q,/uj 
(as opposed to the one-dimensional case) is that the ra- 
dius of the toroidal trap corresponding to V e g increases 
with Q/uj. 

We also note that with increasing at fixed era, cen- 
ter of mass excitation always occurs before quantized 
vortices are formed. The states ^ C m and ^ v cannot be 
connected by low-order application of any single-particle 
operator because of their small overlap, as noted above. 
Thus, even in the region where ty v has a lower energy 
than f cm , the system can remain in ^ cm for times which 
can be very long (i.e., that scale at least linearly with the 
number of atoms N.) The calculation of this time scale 
remains a difficult open question. On the other hand, if 
one first rotates the gas above the condensation temper- 
ature and then cools down, the system will end up in the 
phases shown in Fig. 2. 



VI. SUMMARY 

We have calculated the phase diagram for a rotating 
Bose-Einstein condensate in a quadratic-plus-quartic po- 
tential when the effective interaction between the atoms 
is effectively attractive. For very weak interactions there 
is a phase of (pure) multiply quantized vortices. In this 
phase the energy of the gas in the rotating frame is mini- 
mized and the vortex states are stable against weak per- 
turbations. For somewhat more attractive interactions, 
there is a mixed phase of single and multiple quantiza- 
tion. For even more attractive interactions, the system 
carries its angular momentum via center of mass excita- 
tion. Finally for even stronger attraction, the cloud can- 
not support itself and collapses to a dense atomic state. 

Author GMK wishes to thank Emil Lundh and Ben 
Mottelson for useful discussions. 
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